IPPP/02/40 

DCPT/02/80 

|hep-pli/0207307 



Photon Distribution Amplitudes in QCD 



Patricia Ball\ V.M. Braun^, and N. Kivel^ 



^IPPP, University of Durham, Durham DHl 3LE, UK 

"^Institut fiir Theoretische Physik, Universitdt Regensburg, 
D-93040 Regensburg, Germany 



Abstract: 

We develop a consistent technique for the calculation of real photon emission in hard 
exclusive processes, which is based on the background field formalism and allows a conve- 
nient separation of hard electromagnetic and soft hadronic components of the photon. The 
latter ones are related to matrix-elements of light-cone operators in the electromagnetic 
background field and can be parametrized in terms of photon distribution amplitudes. We 
construct a complete set of photon distribution amplitudes up to and including twist-4, 
for both chirality-conserving and chirality-violating operators. The distribution ampli- 
tudes involve several nonperturbative parameters and, most importantly, the magnetic 
susceptibility of the quark condensate. We review and update previous estimates of the 
susceptibility and also give new estimates of parameters describing higher-twist amplitudes 
from QCD sum rules. 
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1 Introduction 



Hard exclusive processes involving photon emission are attracting increasing attention. 
Examples include transition form factors like 7* 717 with one real and one virtual 
photon deeply- virtual Compton scattering (DVCS) and rare radiative i?-decays 
like B — i> £z/£7, B frf and B K*^ 0. A specific feature of the QCD description of 
such processes is that a real photon contains both a "pointlike", electromagnetic (EM), 
and a "soft", hadronic, component. This distinction is familiar from the studies of the 
photon structure functions in deep inelastic scattering (see e.g. [|]), but so far has not 
been studied, in a systematic way, for exclusive processes. It is the objective of this work 
to develop an efficient formalism for describing either component of the photon and to 
give an update of the available information on the hadronic one. 

We will work in the background (BG) field formalism. This technique was pioneered by 
Schwinger for QED |^ and later rediscovered as a convenient tool in QCD, in particular in 
connection with calculations of the QCD /5-function 0, power-corrections to deep-inelastic 
scattering and QCD sum rules p, ^. An excellent technical review can be found in 
T0| . The principal advantage of the BG field approach is its explicit gauge- invariance 



which allows the use of different gauges for the soft (classical) and hard (quantum) fields 
and also avoids contact terms; instead of the latter, the formalism involves new operators 
containing the EM instead of the gluonic field-strength tensor. It also allows an intuitive 
physics interpretation which is important by itself. As an example , consider quarks and 
antiquarks in the QCD vacuum, in a constant (electro)magnetic field. In the weak-field 
limit, the induced magnetisation of the vacuum is proportional to the applied field, the 
quark density, the quark electric charge Cq and the parameter x, the so-called magnetic 
susceptibility of the quark condensate: 

{Q\qcrci3q\Q)F = eqx{m) Fai3- (1-1) 

Here (gg) is the quark condensate, Fap the field-strength tensor of the external EM field, 
and the index F indicates that the VEV is taken in the vacuum in the presence of the 
field Fafi- If the magnetic field is allowed to vary with a certain frequency, the response 
becomes sensitive to the quark-antiquark separation and gets more complicated. In the 
limit of light-like separations = 0, which is relevant for hard exclusive processes, the 
magnetic susceptibility has to be substituted by a response function 0^ ||TT 

1 

(0|g(z)o-a/3g(-z)|0)F = egx{qq) J duFapUl - 2u)z)(j)^{u) 



1.2) 



where the normalization is chosen such that du(j)^{u) = 1. For a plane-wave configu- 
ration Fai3{z) ~ exp{—iqz), q^ = 0, and in the infinite- momentum frame g+ — > 00, the 
matrix-element on the l.-h.s. of (|1.2| ) describes the probability amplitude that a real photon 
with momentum q dissociates into a quark-antiquark pair at small transverse separation. 
The function ^^(m) can thus be identified with a photon distribution amplitude (DA) with 
u and 1 — u being the momentum fractions carried by the quark and the antiquark in the 
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photon, respectively, in full analogy with DAs of mesons |T2 



photon DAs have been introduced and studied for the first time in Refs. |11 



To the best of our knowledge. 
Some more 



results on higher- twist distributions can be found in |jl3[, while in 0] the leading- twist 
amplitude has been calculated using a chiral quark model in the instanton-vacuum. 

The subject of the present paper is the detailed study of these DAs, up to and in- 
cluding twist-4. Most of the discussion will refer to asymptotic real photon states, which 
correspond to a plane-wave configuration of the EM BG field, whose frequency is, however, 
not assumed to be small. This deviates from the usual procedure of expanding in slowly 
varying (classical) fields and requires special techniques which have been worked out in 
Ref. 

The paper is organized as follows: Section 2 is mainly introductory and summarizes 
notations and basic features of the BG field method. Section 3 contains a detailed discus- 
sion of the leading-twist photon DA. Section 4 is devoted to photon DAs of twist-3 and 4. 
Finally, in Sec. 5 we summarize. The paper also contains several appendices with results 
of more technical nature. 



2 The Background Field Method 
2.1 Gauge- invar iance and contact terms 

The basic idea of the BG field method is to modify the quark part of the QCD action by 
including an EM field -B^ in addition to the colour gluon field A^. In order to preserve 
EM gauge-invariance, the covariant derivative becomes 

D^ = d,~tgA;^t''-tegB^, (2.1) 

where = eQg = y/Ana Qq with Qu,c,t = +2/3, Qd,s,b = ~l/3 is the quark electric charge 
and a = 1/137. . . is the fine-structure constant. The action reads (for one flavour) 

S, = jd^xq{x)[iP -m,]g(x) = SQC° + e,y"d^x5^(x)j^(x), (2.2) 

= qjfiq is the vector current. The EM field is treated as purely classical, i.e. it 
does not participate in loops, and is weak, i.e. we only consider expressions linear in i?^. 
In Lorentz- gauge, an ingoing photon with momentum q corresponds to the plane-wave 
field-configuration 

5,(x) = eWe-^''^ (2.3) 

Here e^^^ is the photon polarization vector with e^^\^ = 0. The corresponding gauge- 
invariant EM field-strength tensor is 

F,,(x)=^(eWg.-g,eW)e-^'?^ (2.4) 

To avoid confusion, we will always write F^jy for the EM and G^u = G^^.t'^ for the gluon 
field-strength tensors, respectively. 

(Almost) any calculation in QCD relies on the separation of scales: the contribu- 
tions of quark and gluon exchanges with large momenta are calculated perturbatively and 
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Figure 1: Schematic representation of a hard-soft separation. 

form a hard contribution (coefficient function), the contributions with small momenta are 
parametrized in terms of various distribution functions, alias operator matrix-elements. 
Such a structure is indicated schematically in Fig. |I|. A real photon (not shown) can be 
emitted from both the hard and the soft block. The advantage of the BG field method 
relies on the fact that it separates hard and soft photon-emission in an explicitly gauge- 
invariant way. This is not necessary, but convenient, as we will see below. 

Technically, the gauge-invariance is seen most directly if, in the calculation of the hard 
block, one uses the following expression for the quark propagator in the BG field^ ||1 5|| : 



S{x) 



i ^ 
2^ 



[x,0]j 



du [x, ux]f{u ^(Ta(i + u(Tap ;?f}egF"^('ux) [mx, 0]i? 



-|- terms containing G, 



a/3; 



(2.5) 



where = is the momentum operator, u = \ — u and 

[x, y\p = Pexp dt {x - y)^,[gA>'{tx + ty) + e,5^(te + ty)] | (2.6) 

is the path-ordered gauge-link. Here and below the subscript F serves to indicate that the 
BG field is included, which means in particular that both the abelian EM and the non- 
abelian colour-phases are present. We will also use the notation [x,y] without a subscript 
if the EM field is not included. 

The terms in Fap generate U{1) gauge-invariant contributions describing photon-emis- 
sion from small distances, whereas the gauge-factors [., .]f in the first term on the r.-h.s. 
of ( |2.5|) eventually combine with quark fields that connect the hard and the soft block, 
producing nonlocal gauge-invariant operators. In the simplest situation with no external 
hadrons, one is left with the vacuum expectation values of such operators in the EM BG 
field, generically 



{0\q{x)r[x,y]pq{y)\0)F = ie,jd^zB^{z) (0|T{g(x)r[x, y]g(y)j^(^)}|0) 





^ For simplicity we present the result for the massless quark propagator. For massive quarks, see 
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+ ze, dt{x - y)^B>^itx + iy){0\qix)T[x,y]qiy)\0). (2.7) 



Here T is an arbitrary Dirac-structure. The first term on the r.-h.s. of ( p.7|) comes from 
the expansion of the action e*^'', cf. (p.2|), whereas the second term originates from the 
gauge-hnk ( [2.61 ). The expression ( |2.7| ) is invariant under U{1) gauge-transformations -B^j 
B^+d^(f). As an example, let us check this for plane-waves. Using (|2.3|) , the r.-h.s. becomes 



ze, / rfSeWe-^'^^(0|T{g(x)r[x,y]g(y)j^(z)}|0) 



+ - ^"'1 (0|g»r[x,y]g(y)|0). (2.8) 

Note that the expression in the first line is nothing but the amplitude of the photon-to- 
vacuum transition 



(0|g(x)r[x,y]g(y)|7(')(g))=^e,yd^zeWe-^«^(0|T{g(x)r[x,i/]g(i/)/(^)}|0). (2.9) 
To check gauge-invariance, replace e|t^'' — >■ g^. The first line in ( |2.8|) then becomes 

e, / d'ze-'^^0\d,T{qix)T[x,y]qiy)riz)}\0), (2.10) 



which is generally nonzero because of contact terms. To calculate the latter, use that 
inside the QCD functional integral the divergence of the vector current can be replaced 
by a functional derivative over the quark (antiquark) field: 

d,nz) = ~q{z) ij^f^') - {j^f^') ^^'^ ■ (2.11) 
Integrating by parts over the quark fields in the functional integral, one obtains 

e, / d'ze-^'^^ {(5(^)(x - z){0\q{z)T[x,y]q{y)\0) - 6(^\y - z){0\q{x)T[x,y]q{z)\0)} 

= e, [e-"^^ - e-^'iy] {0\q{x)T[x,y]q{y)\0) . (2.12) 

This is exactly what is obtained from the second term in ( p.8| ), but with the opposite sign. 
One sees that restoration of U{1) gauge-invariance is achieved through elimination of the 
contact terms. 

Calculation of the hard scattering amplitudes (hard block in Fig. |l|) can usually be 
done for on-shell partons. In the operator product expansion language, this means that 
operators which vanish by virtue of the QCD equations of motion (EOM) can be omit- 
ted. The case of real photon emission is special, because the EOM operators can generate 
nonzero contact terms in correlation functions with the EM current. The BG field tech- 
nique provides an important simplification, since no contact terms are present whatsoever. 
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but instead the EOM are modified by the presence of the BG field. As a simple example, 
consider the following local operator with two covariant derivatives: 

Of,y = g7M7a DaDu q ■ (2-13) 
For hadronic matrix-elements, one can use the identity [Da,D^] = —igGau to write 

{h'\0^,\h) = -i{h'\q^^^^gG^,q\h) + {h'\q^^ Du^ q\h) (2.14) 

and neglect the last term. For photons, one either has to take into account the EOM 
operator explicitly, or use the BG field formalism, in which case there is an additional 
term in the commutator: [Da, Dy\ = —igGau — ieFau- We get 

(0|C^.|0)i. = -i{0\qj^jagGauq\0)F {qq) + {0\qi,DuP q\)F. (2.15) 

The third term on the r.-h.s. is again zero by virtue of the EOM, p q = 0, but there is 
an additional contribution containing the photon field (and the quark condensate {qq)), 
which illustrates the generic feature of the BG field formalism, the appearance of new 
operators containing the EM field-strength tensor. 



2.2 Light-cone expansion 

After these general remarks, we proceed to the explicit construction of the light-cone 
expansion of quark- ant iquark operators in the BG field. We assume that the quark and 
the antiquark are separated by the distance and take the limit ^ keeping the scalar 
product {q- x), q"^ = 0, fixed. As discussed in detail in [0, |18[, the light-cone expansion of 



a generic nonlocal operator produces two sequences of terms that are singular and analytic 
on the light-cone, respectively:^ 

{Q\q{x)T[x,y]Fq{ymF = (0|g>)r[a;, 2/]^g(i/)|0)f ^'^^^'^ + (0|g(a;)r[x, y]pg(y)|0)r''*^ 

(2.16) 

The singular contributions correspond to small distances, i.e. to the coefficient function 
in the OPE expansion, and the analytic contributions induce the corresponding operator 
matrix-elements. The separation between singular and analytic contrubutions involves, as 
usual, a factorization scale ^p- In the case at hand, the singular contributions are given 

by 

mx)T[x,yU{y)\Q)T''''^' = -Tr [T[x,y]FS{y - x)] , (2.17) 

where the trace is taken over both spinor and colour indices and S{y — x) is the quark prop- 
agator in the BG field ( p.5| ). As for the analytic contributions, they can be parametrized 
in terms of nonperturbative photon DAs of increasing twist. To leading order in the QCD 
coupling and for massless quarks, one obtains for the chiral-even operators 

1 

{0\q{x)jf.[x, -x]Fq{-x)\0)F = jdu x" e,F,^{-ix) {-^^ ^ + G^'^H^) + 0{x^) 



'-Gi. in particular Eqs. (3.7)-(3.11) in M 
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(0|g(x)7M75[a;, -x]Fqi-x)\0)i 



<x) 



Nr 



+ G^''\u) + 0{x^)\i2.l%) 



whereQ 

and G^'"\u)^ G^"'\u) are nonperturbative functions that depend on at most logarith- 
mically and will be determined later. Here and below we use the standard notation 
[1 /2)eaf3fiuF'^'^ and the conventions for the e tensor and 75 matrix as defined in |jl9|. 



a/3 



Counting dimensions, one finds that the perturbative contribution in (|2.18|) is twist-1, and 
the nonperturbative corrections are twist-3. Note that this dimensional twist-counting 
refers to the "dynamical" twist of a matrix-element, as opposed to the "geometric" twist 
of a (local) operator. We will discuss this in more detail in Sec. 4.1. 

For chiral-odd operators the perturbative (singular) contribution vanishes identically 
for massless quarks. One then has |T3[ 



\q{x)af,^[x,-x]Fq{-x)\0)F = eg{qq) J duF^y{-ix) 



+^eq{qq) j dux^ {x^F^p - x^F^p} (-^ 


{0\q{x)[x, -x]Fq{-x)\0)F = 0. 



X B M + ... 



(2.19) 



The function 0^(-u,/i) is the (nonperturbative) photon DA of leading twist-2, while the 
terms in A{u) and B(m) correspond to twist-4 corrections. They will be discussed in detail 
in the following sections. Note that the vacuum expectation value of the scalar operator 
does not contain any contribution linear in the BG field and thus vanishes identically. 

For the strange quark, and also considering possible lattice calculations of DAs, it 
can be interesting to keep the lowest order term in the quark mass m^. To leading- twist 
accuracy one then has 



{0\q{x)af,u[x, -x]Fqi-x)\0)F 



J duegFp^{~^x)l - ^"^"^ 



47r2 



ln[x^/i^M(l — u)] + const 



(2.20) 



(cf. App. A in |jrT[). Note that in this case there is a weaker, logarithmic singularity on 
the light-cone. The calculation of the constant under the logarithm requires a careful 



matching with the UV-divergent contribution of the quark loop to the photon DA |11 



■^The definition of the integration variable u is a convention. It is chosen such that in the case of an 
ingoing photon u corresponds to the momentum fraction carried by the quark. 
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It is instructive to verify that the perturbative (singular) contributions to ( p.l8| ) indeed 
correspond to the perturbative wave-function of the real photon in the infinite-momentum 
frame. To this end choose plane- wave, (|2.4| ), and the coordinate system such 

thatj^ q = {q'^,0,0), x = (0,x_,r), x"^ = — r^, where r is a two-dimensional vector in the 
transverse plane. For the leading- twist "plus" component, assuming physical (transverse) 
photon polarization, one obtains 



(0|g(x)7+^[x, -x]M-mF = --^.q+ldue'^'^^^' 



e(^) •r)(2n- l)±ieifcrie 



(A) 



(2.21) 

where e^fc = eik+~ is the two-dimensional antisymmetric tensor in the transverse plane, 
i,k = 1, 2, which is the expected result. Indeed, the photon wave-function can be defined 

as 



i^hiu, r) 



dx- 



71 



iNcCq 



i2;-(j+(2u— 1) 



(0|g(x)7_ 



X, -x]Fq{-x)\0)F 



x+=0,x I =r 



|(eW-r)(2?z-l) + /ize,fc e^^^} 



(2.22) 



where h = ±1 corresponds to the quark helicity. Choosing the photon polarisation vector 
(0, 1, ±i, 0)/v^, one obtains the familiar expressions given in, for instance. 



as e 



(±) 



©HQ. 



3 The Leading-Twist Photon Distribution Amplitude 

We proceed with the detailed study of the leading-twist DA of the photon. Like for a trans- 
versely polarized vector meson p2|, the leading-twist contribution corresponds to the 
chiral-odd Lorentz- structure. Following ||Tl| we define the DA as the vacuum expectation 
value (VEV) of the nonlocal quark- ant iquark operator with light-like separations, 

1 

{0\q{z)aap[z, -z]Fq{-z)\0) F = egx{qq) / du Fc,p{-^z)(f)^{u) , (3.1) 



where z"^ = and we have inserted the path-ordered gauge-factor which was implied, but 
not shown in ( [1.2| ). This definition is equivalent to 

1 

t jd'ye-"'y{0\T{f{y)q{z)a^p[z,-z]q{-z)}\0) = i{qpg - q^Q ,p)x {m) j due'^'^'^.iu) 



+ ^ [e''' - e-*^^] {0\qiz)a^p[z, -z]qi-z)\0), (3.2) 
qz 



"'We define = (a*^ ± a^)/^/2 for any four- vector a. 
^Up to a different normalization. 
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cf. ( p.7| ) and ( p.8| ). The first term on the r.-h.s. of ( |3.2| ) corresponds to the VEV in the BG 
field of the plane-wave, ( |3.1| ), ( p.4| ), and the second term is the contact term. In fact, in the 
above case the contact term vanishes: the VEV {0\q{z)aaf3[z, —z]q{—z)\0) does not include 
any external momentum and from Za alone it is not possible to build an antisymmetric 
Lorentz-structure (in (a,/3)). 

Yet another form of the definition ( p.2| ) is (cf. ( |2.9| )) 



{0\q{z)ac,f3[z,-z]q{-z)\-f^^\q)) =tegx{qq) (qpel^^ 



(3.3) 



Here we assume that the physical photon polarization is transverse to the (g, z) plane. 
The three definitions (|3.1| ), ( p.2|) and are equivalent and can be useful in different 
contexts. 

As always in quantum field theory, the extraction of the asymptotic behaviour — > 
creates divergences that have to be regularized. As a result, both the magnetic suscepti- 
bility X aiid the DA 4>^{u) are scale-dependent. To leading-logarithmic accuracy: 



6uu 



1+ Yl L^"^"'PrMCT{u-u) 
n=2A,... 



where Cn^'^{x) are Gegenbauer-polynomials, b = ll/3Nc — 2/3nj, L 
jgq = —3Cf is the anomalous dimension of the quark condensate and 

n+l 

7n = Cf 



n+l 



(3.4) 

as(At)/as(Ato), 
(3.5) 



are the anomalous dimensions of the chiral-odd local operators of leading- twist [Q. 

Photon DAs are interesting mainly because of the fact that the magnetic susceptibility 
X appears to be rather large. A crude estimate of its value can be obtained in the vector- 



dominance approximation [25|: consider the correlation function 



I / Ae-"'^(0|T{j^(y)g(0)a„;3g(0)}|0) = i{qpg,^ - qag,p)x{q', f^) (qq) 



(3.6) 



where in contrast to ( p.2|) the separation between the quarks is set zero and can be 
arbitrary, instead of = for a real photon. From the normalization of the photon DA 
it follows that x(/^) = x(?^ = 0,/x). Assuming that the correlation function is dominated 
by the contribution of the lowest-lying states po and u and using the standard definitions 
of couplingsp] in the flavour-SU(2) limit 



{Om,qY='^'^\poico))=ei'^fvmv 



215 MeV [26], 

(160± 10)MeV |2§ , (3.7) 

®The currents with specific isospin are given by [qTqY^^^^"^ — (uFu =F dTd)/\/^. my is the average 
mass of Lu and p(770). 



qa)f^ 



V ' 



fv 
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one obtains 

(3.8) 

my — q'' 

and therefore ^ 

XVDM ^ -- - 2.7 GeV-2, (3.9) 
mv{qq) 

using the value (qq) = — (250MeV)'^ (at 1 GeV)[]. Ahernatively, one can assume that the 
vector- dominance approximation ( |3.8|) is vahd for large Euclidian q^ —>■ — oo, in which limit 
the leading contribution to the correlation function (|3.6| ) is given by the quark condensate: 

X(g^) =-°° A ■ (3-10) 

-r 



Comparing this expression with ( ^.81 ), one obtains an even simpler estimate 

2 



which is rather close to ( p.9| ). Here the subscript "LD" stands for the so-called local-duality 
approximation. 

A more sophisticated study including contributions of higher resonances to (|3.8|) on one 



side, and contributions of higher- dimensional operators to (|3.11|) on the other side, com- 
bined with the usual QCD sum rule techniques to match between the two representations 
in the most efficient way, was done in E7\, 53. The result of this analysis. 



XsR(lGeV)~4.4GeV-^ (3.12) 

was used in numerous QCD sum rule calculations of EM hadron properties over the last 
15 years. In App. B we give an update of the analysis of using new data for the 

resonances, the most recent value of as and including 0{as) radiative corrections to the 
sum rule. We obtain 

Xsr(1 GeV) ~ (3.15 ± 0.3) GeV'^ , 

/^^(l GeV) = xm = -(50 ± 15) MeV, (3.13) 

which is smaller compared to the older estimate, mostly due to the large negative radiative 
correction. 

The contribution of the photon DA (twist-2) to a generic hard exclusive process is, 
compared to the perturbative contribution (twist-1), suppressed by one power of the hard 
scale Q and is of order 

1 (27r)^|(gg)| 650 MeV 

^ Q - ^ ' 

This relatively large mass scale is due to the large magnetic susceptibility and to the well- 
known fact the the quark condensate contribution to physical observables is enhanced by 



^The sign of x depends on the sign convention used for the EM couphng in (2.1). 
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a factor (27r)^ compared to the perturbative contributions because of a different structure 
of the phase-space integrals, cf. Eqs. ( |2.18|) and ( |2.19| ). 

The available information on the shape of the photon DA is more controversial. At 
large virtualities, the shape is uniquely determined by the renormalization properties of 
the relevant operators, cf. ( |3.4|) . As the anomalous dimensions in (|3.5|) rise with n, one 
obtains, for very large scales fj,, the asymptotic photon DA 



/.asy 



[u) = 6u{l-u), (3.15) 



which is the same as for it, p, . . . mesons. In Ref. arguments were given suggesting 



that at smaller virtualities of order 1 GeV the photon DA is still close to the asymptotic 
form. The physical picture behind the argumentation in |jTl| was that corrections to the 



asymptotic DA are nonperturbative in nature, and are expected to be washed out when 
one sums over contributions of different hadron resonances p, p', p" , .... This qualitative 
conclusion was supported by the calculation of the second Gegenbauer-coefficient 02 at the 
scale 1 GeV using QCD sum rules, yielding a small number. Unfortunately, this estimate 
appears to be very sensitive to the choice of input parameters, see App. B. In addition, 
there has been increasing evidence in the last years indicating that "standard" QCD sum 
rule estimates of corrections to asymptotic DAs are not reliable. We therefore believe that 
the existing evidence in favour of the asymptotic DA of the photon at low scales is only 
qualitative. 

Another calculation of the photon DA was done in using the instanton- model of 



the QCD vacuum. This approach neglects altogether perturbative gluon radiation which 
generates the asymptotic DA, and the results correspond to a very low normalization point 
of order of the characteristic instanton-size pj ~ 600 MeV. The resulting photon DA can 
to a good accuracy be represented by a simple expression 



M, = 600 MeV) = [a-6u{l-u)] (3.16) 



1 



with a ~ 7.58.0 In the instanton-model the shape of the pion and the photon DA turn 
out to be very different. The reason for this is that for the pion there exists an additional 
(to perturbative radiation) mechanism that suppresses the DA at the end-points P |. 



Roughly speaking, this effect is due to a nonvanishing size of a constituent quark when 
it is probed by the axial current, which is required quite generally by chiral symmetry 
breaking. For the vector current (photon) the constituent quarks appear to be pointlike, 
and the suppression is not required. As a result, the pion distribution is predicted in this 
model to be close to the asymptotic one already at very low scales, while the photon DA 
is almost flat and does not vanish at the end-points, see Eq. ( ^.16 ) and Fig. 0. From our 



point of view this evidence is not conclusive either, since the observed end-point behaviour 
corresponds to contributions of large invariant masses, and it is not clear where inclusion 
of them in the instanton-model calculation is legitimate. By its physical meaning, the 
distribution amplitude corresponds to the contributions of low momenta and we expect 
that an asymptotic-like shape would be recovered if the calculations were made with an 



^V.M.B. thanks M. Polyakov and C. Weiss for providing this parametrization. 
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Figure 2: The chiral-odd photon DA of leading twist. The two curves correspond to 
the asymptotic DA ( |3.15| ) (sohd) and the instanton model at the low normalization point 
0.6 GeV ( CT ) (dashed). 



explicit ultraviolet cutoff of order of the instanton-size. We believe that this question will 
eventually be decided by lattice calculations and/or experiments on hard exclusive dijet 
photoproduction in which the shape of the photon DA can be measured ||29|| . 

For the strange quark, one might want to take into account the quark mass correc- 
tions, in which case there is the complication that the separation of singular and analytic 



contributions on the light-cone involves a factorization scale, cf. ( 2.20 ). For completeness 



we quote here the corresponding result from [11 



{0\s{x)aap[x, -x]s{-x)\'y^^\q)) = ie^ 




27ms 

2 

9^ 



(u - 



3ms 
4^ 



w 2 -^ 367r2 ^ 
in[—x WqUU) 5 — h 2'yE 

4 



(3.17) 



The quantity Wq = 2 GeV^ is the continuum threshold in the classical QCD sum rule 
for the 0-meson [0] and An'^/g'^ ~ 11.7 is the 0-meson coupling (squared) to the EM 
current. By definition, the magnetic susceptibility of the quark condensate does not include 
perturbative contributions of the quark loop ~ m^. 



4 Distribution Amplitudes of Higher Twist 
4.1 General classification 

The aim of this section is to give a general classification of photon DAs up to twist-4. As 
already stated in Sec. 1, we define DAs as (renormalized) matrix-elements, in the EM BG 
field, of nonlocal gauge-invariant operators at strictly light-like separations. For brevity, 
starting from this section we will not show path-ordered gauge-factors between the fields, 
they are always implied. We distinguish between matrix-elements of operators with an 
odd number of Dirac 7-matrices, which we term chiral-even, and operators with an even 
number of 7-matrices, which we refer to as chiral-odd. 

The twist-counting employed in this classification is based on counting of powers of 



the photon momentum in the infinite- momentum frame, cf. the discussion in To 
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T wi st 




( iiu) OCT n 


2 




+ -L (f)j 


3 






4 




± — 



Table 1: Identification of two-particle DAs witli projections onto different liglit-cone 
components of tlie nonlocal operators. For example, _L refers to ^7^75(3'. 



Twist 








3 


+ ±+ A V 






4 




Ti 

+ ±±± T2 
+ - + ± T3 
+ T4 


+ ± 5 5 



Table 2: Same as Tab. [1| for three-particle DAs. -|- _L_L refers to qG^±'j±j5q etc. 

explain this, take a plane-wave configuration for F^^, and choose the frame of reference such 
that q^ = (g+,0, Oj_), 2;^ = (0,2;_,0_l) with g+ 00. The twist-expansion corresponds 
to counting powers of the large momentum g+, taking into account that the photon's 
transverse polarization vector is of order e^i^'' ~ {q+)^ and ~ Matrix-elememts 
of nonlocal operators typically involve several Lorentz-structures which can be separated 
by suitable projection operators. In this way, the leading-twist structures correspond to 
the maximum number of "plus" components. Each replacement of a "plus" projection by 
a transverse one increases the twist by one unit, and each "minus" projection adds two 
units of twist. 

The complete classification of the relevant two- and three-particle DAs and their re- 
lation to specific light-cone projections of the matrix-elements is given in Tabs. 1 and 
2, respectively. As is well known, "minus" components of quark-field operators do not 
describe independent partonic degrees of freedom, but can be expressed in terms of mul- 
tiparton DAs. The explicit definitions and relations between various distributions will be 
worked out in the next sections. 

Note that the classification of photon DAs is very similar to the classification of the DAs 



of the transversely polarized p-meson discussed in |^T], For the twist-3 distributions 
considered in the next section the analogy is in fact one-to-one. For twist-4 distributions 
the situation is somewhat more complicated because of the different structure of the QCD 
EOM in the presence of the BG field. 
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4.2 Twist-3 distribution amplitudes 

The twist-3 DAs of a real photon are related to matrix-elements of chiral-even operators. 
Since the vacuum expectation value of the vector current in the BG EM field vanishes 
by the EOM, (0157^^10)^ = 0, the relevant local operator with lowest dimension has to 
include the gluon field: 

{^\qgG^uiai^q\^)F = e^h^Do^F^^, (4.1) 

where Da is the covariant derivative. The nonperturbative constant /s^ has dimension 
GeV^ and provides a natural mass-scale for the twist-3 DAs. There are two two-particle 
DAs which we define as (cf. Tab. 1) 

(0|g(^)7/.g(-^)|0)F = ej37 I dmP^''\u,fi)F,,{-^z), 



-1 

{0\qizh>.l6qi-zmF = e,h,'- I du^^''\u,fi)F,^i-^z) , (4.2) 







where 



F = F 



We will use similar shorthand notations in what follows. The normalization of the DAs 
is such that du'ijj^'"\u) = duip^"'\u) = 0. Note that the regularization (and renor- 
malization) removes the singular l/z"^ contributions, so that the matrix-elements in ( [4.2| ) 
correspond to the analytic parts of the amplitudes in the light-cone limit, in the sense of 
the separation in Eq. ( |2.16|) . Comparing (|4.2| ) with ( p.l8| ) we find 



G(^)(«) = /3,V^(^)(w,/x~l/|x|), 

G^^\u) = ^h,^'~^\u,fi^l/\x\). (4.3) 
In addition, there exist two three-particle chiral-even twist-3 distributions: 

{0\q{z)gG^^{vz)'ja'j5qi-z)\0)F = Cqfs^ j VaA{a)DaFf,y{ayz), 

{{)\q{z)gG^^{vz)i-faq{-z)\^)F = ejs^ j VaV{a)DaF^^{a^z) , (4.4) 

where ay = aq — Uq + vOg, and a is the set of the three momentum fractions a = 
{aq }. The integration measure is defined as 



Va = I daq / dag / dag 5(1 — ) ■ (4-5) 

Jo Jo 



Alternatively, the same distributions can be defined as matrix-elements between the ingo- 
ing photon state and the vacuum, and in this form they are very similar to the definitions 
adopted in [^l[ ^ for the transversely polarized p-meson, apart from a different overall 
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normalization. For two-particle distributions: 



1 



\q{zh,l,q{-zW>{q)) = -eJ,,e,,,Jll du e'^^^ ij^^\u, fx) . (4.6) 

The functions ip^^^u.fi) and ip^'^^u, fx) are related by integration by parts: 

= 2 / . (4.7) 

Jo 

For the three-particle distributions one has 

(0|g(;.)(7G^,(t;z)7„75g(-^)l7^'ng)) =ej37?a[9.e^^-g^e^i] J VaA{a)e-'^^'^% 

(0|g(^)(7G^,(t;^)z7ag(-;^)l7^'H9)) = ej37ga[g.e?^-g^e2] / PaV(a)e-^^^"^ (4.8) 



Like the corresponding p-meson DAs, the two-particle distributions ip^^\u) and il)^°^{u) are 
not independent and can be related to three-particle distributions using the QCD EOM. A 
closer inspection shows that all corrections to the EOM in the BG field vanish in the chiral- 
even case, so that these relations remain intact. | As result, the photon DAs can be read 



off the corresponding expressions for the p-meson ||3^ , where one must omit the so-called 
Wandzura-Wilczek parts that are due to the nonzero meson mass. To next-to-leading 
order in the conformal expansion, cf. Sec. 4.3, we obtain 



V(a) = 540a;}f(aq — aq)aqaqa 



A{a) = 360 aqa^a^ 



9 ' 



1 + - (7a, - 3) 



5(3^' - 1) + ;^(15cu^ - 50(3 - 30e'-^ + 35^') 



(1 -a(5e^- 1)^(1 

-20^i< + }^{uj^- 3.;^) uu$,{7e - 3), 



(4.9) 



where the parameters oj^ and u]^ correspond to matrix-elements of local operators of 



dimension 6, see ||3^ and App. B. 

In vector-dominance approximation one finds 



VDM 

37 



-/;C3^-(4±2)-10-^GeV 



, ,y,VDM 
'^7 



UJ. 



,A,VDM 

7 



V 



3.8 ±1.8, 
-2.1 ± 1.0 



(4.10) 



^This is in contrast to the chiral-odd ones, to be discussed in the foUowing subsection, where one VEV 
is nontrivial, and induces a nonzero conformal expansion of an exactly vanishing DA. 
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where /pCs — (20 ± 10) MeV (see App. B) is the p-meson couphng to the quark-antiquark- 
gluon current in ( [4.1|) . Our result for this couphng is by a factor 3 larger than an old 
estimate in obtained in Ref. [Q, and quoted in |^T], |32|. We have estimated corrections 
to the vector-dominance approximation for /s^ using QCD sum rules, see App. B, with 
the result that these corrections are smaller than the (large) error bars for the p-meson 
contribution. We, therefore, accept the vector-dominance approximation for the present. 
The scale-dependence of the twist-3 parameters is given by, cf. Ref. pT], {Ca = Nc): 



and 



/3>2)=L^'/%(/i^), y = -ic^ + 3C^, (4.11) 



r>r^ 2/^ 2/^ 2/^ 

5/^ iri Ir-i I r-i 

— oOyl nO_p -|- 



(4.12) 



4.3 Twist-4 distribution amplitudes 
4.3.1 Definitions 

Twist-4 DAs, which are all chiral-odd, are more involved and require a detailed discussion. 
The general parametrization of the nonlocal quark-antiquark operator on the light-cone 
involves two invariant functions: 

1 

{Q\q{z)af,yq{-z)\{])F = eqx{qq) j du(j)^{u, fx) F^^{-^z) 



1 

+ ^eg(gg) J du h-y{u, /i) {z^Ff.^ - z^F^^} (-^z) , 



{0\q{z)q{-z)\0)F = 0. (4.13) 

Here (p^yiu) is the leading-twist distribution discussed in Sec. 3 and h^{u, fi) is a new DA 
of twist-4, cf. Tab. 1. Note that the third possible Lorentz-structure z^ {D^F^p — D^F^p} 
is not independent and can be eliminated using the QED Bianci-identity. 
The three-particle DAs are more numerous and can be defined as: 



{0\q{z)gG^^{vz)q{-z)\0)F = eg{qq) J VaS{a)F^^{a^z) , 
{Q\q{z)gG^,{vz)i^^q{-z)\())F = eg{qq) J VaS{a)F^,{a,z) (4.14) 



and 

{Q\q{z)aapgG^y{vz)q{-z)\{])F 
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+ iegiqq) [gap{gp^D^ - g^^D^} - {a ^ {3)] j Vaf2{a)Fp, 
+ ^eg(gg) [Df^z^ - D^z^] j VaT^{a)Fap{ayz) 

+ i(iq{.m) [DaZfi - DfsZa] j 'Dafi{q)F^y{ayz) , (4.15) 

see Tab. 2 for the identification of the relevant hght-cone projections. 

In addition, we have to consider also quark-antiquark-photon operators: 

{Q\q{z)e,Fp,{vz)q{-zmF = e,{qq)F^,{vz) (4.16) 

and 

{0\q{z)a^pegFp,{vz)q{-z)\0)F = 0. (4.17) 
It is convenient to rewrite these expressions introducing auxilary three-particle DAs: 

{0\q{z)e,Fp,{vz)q{-z)\0)F = e^iqq) jvaS^{a)F^,{a,z), (4.18) 

{0\q{z)aaf3egFf,^{vz)q{-z)\0)F = ieq{qq) [D^zp - DpZo\ j Vaf2{a)Fp^{a^z) + . .(.4.19) 
with 

S^{a) = 6{aq)6{aq), f2{a) = . (4.20) 

As explained in the next subsection, we need the conformal expansion of S-y{a) and T2(a), 
which turn out to have nonzero coefficients. 

In some applications it can be advantageous to rewrite the definitions ( [4.13| ) to ( [4.19| ) 



for the particular case of a photon in the initial state. For two-particle DAs we obtain 

1 

{0\q{z)a^pq{-z)\^^'\q)) = le^x^qq) (g/^e^ - qJ^A [ due'^^^ <f>,{u, f,) 



1 . {qq) 





1 



+ -teq ^ [zpe^^^ - zjp^') I du e^^"^ h^{u, fx) . (4.21) 







It is easy to see that 

h^{u) = —A I da{u — a)h^{a, fi). (4.22) 







C-parity implies the symmetry h^{u) = —h^{l — u) and h^{u) = h^{l — u) 
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For three-particle distributions the corresponding definitions read 



and 

m{zW^9G,,{vz)q{-z)\^'^^\q)) = 

e.(gg)[g«ei%^. - g,e?^,i - g^e?^,,^, + g,eK,i]T,(., ,.) 

+ ^[g.g.ei>. - q,q,e^^lz. - q^qj^z, + q,qJllz,]Uv, qz) 
qz 

+ ^[g.g.eK., - g,g,e?).„ - g.g.e?>, + g,g.e?>jT,(., g.), (4.24) 



g2; 



where 



Ti{v,qz) = j Vae-'^'''^"-Ti{a). (4.25) 

The functions Tj(a) are related to the distributions for a generic EM BG field by integration 
by parts: 

Ti{a) = —2 / duTi{u, aq + aq — u) . (4.26) 
Jo 

Here we imply Tj(a) = Tj(ag, a^). The corresponding equivalent definitions of S*^ and Tj 
are analogous to ( [4.23| ) and ( |4.24| ), respectively; the relation between Tj and Tj is the 



same as in ( ^.26| ). 



4.3.2 Conformal expansion of three-particle twist-4 distribution amplitudes 

In practical calculations one has to use models for the higher-twist DAs with a few non- 
perturbative parameters. The major difficulty in constructing such models is that the 
QCD EOM imply the existence of relations between different amplitudes and that these 
relations have to be obeyed identically. The standard approach |jll], |3^ is to expand the 
amplitudes in contributions of operators with increasing conformal spin and retain only 
a few leading terms. The reason why this helps is that the QCD EOM are essentially 
tree- level relations that respect the same symmetries as the QCD Lagrangian. Relations 
between DAs imposed by EOM are therefore "horizontal" in the sense that they relate 
contributions of only the same conformal spin. This implies that the EOM can be solved 
order by order in the conformal expansion and, most importanly, the truncation of the 
conformal expansion at a certain order is not in conflict with the QCD EOM. In addition, 
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as familiar from leading- twist, the conformal expansion is consistent with renormalization 
at LO level: contributions with different conformal spin do not mix with each other to 
leading-logarithmic accuracy. This ensures that higher-spin contributions omitted in the 
model at a certain scale will not reappear at different scales through the evolution. 

Technically speaking, the conformal expansion of distribution amplitudes is the ex- 
pansion in orthogonal polynomials that correspond to irreducible representations of the 
conformal group. A suitable basis is defined in App. A. We have already seen examples 
of conformal expansion in Eqs. (|3.4]) and ( |4.9| ). The expansion of the three-particle DAs 
of the photon of twist-4 is completely analogous to that of the p, done in ||32| . Without 
going into details we just quote the results: 



Ti{a 
^2 (a 
^3 (a 

S{a) 



-120(3(2 + C2^)("g- - ag)agagag, 

30al{ag - ag) - + (Ci - O (1 - 2«3) + (2(8 - 4«, 
-120(3(2 - C2^)(ag - ag)agagag, 

SOaliag - ag) [{k + k+) + (Ci + (t) (1 " 2«<;) + (2(8 - 4a, 
30«J{(fi: + (1 - ag) + (Ci + Ci')(l - ag){l - lag) 



+C2[3(ag - 
-30aJ|(K 

+C2[3(ag - 



ag{\ 



0L„\^ - 



Otr. 



ag{l 



"9)]} ' 

h(Ci-0(l-a3)(l-2aJ 
)]}• 



a 



(4.27) 



Here terms in k and correspond to the lowest conformal spin-3 and the contributions 
of Cs are next-to-leading with conformal spin-4. For the lowest spin the definitions are: 



{Q\qgG ^yq\Q) F = eq{qq){K^ + k) F^^, 
{^\qgG^,uim\^)F = eq{qq) (k+ - k) F^^; 
K and renormalize multiplicatively [p!l|] : 



(4.28) 



«:+(p2) = L(^^->-^)/''«:+(p^), = 2>Ca 
/^(/i^) = L(^-->-'/^)«(/i2), 7_ = AC A - 



— — Cp, 
3 ' 

3Cf. 



(4.29) 



The ^-parameters can be expressed in terms of reduced matrix elements ((Q'-*-')) of lo- 
cal quark- ant iquark-gluon operators with one extra covariant derivative and also receive 
contributions of quark-antiquark-photon operators. One finds 



Ci 



^2 = l + l((g(3)))_^ 
22 22^^ 220 

7 



{{Q 



(4.30) 
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Ci 


C2 






0.2 





0.4 


0.3 









Table 3: Numerical values of parameters of twist-4 chiral-odd DAs taken from [11]. Renor- 
malization scale is n = 1 GeV. 



where Q^^\ Q^^\ and Q^^^ are multiplicatively renormalizable operators (in pure QCD, 
that is not including mixing with photon operators) listed in Sec. 4 of Ref. |]Tl| : 



with 7q(i) 



2 ^ 



7q(5) 



5C, 



(4.31) 



The exact definition of these operators is not important here. Numerical estimates ||Tl 
are collected in Tab. |^. 

Note that four parameters are defined by three independent matrix elements so that 
there is one (exact) relation: 



Ci + 11C2 - 2C^ 



7 
2 



(4.32) 



which is a consequence of the theorem by Ferrara-Grillo-Parisi-Gatto |35|, which states 
that the divergence of a conformal operator vanishes (in a conformally- invariant theory). In 
the present case, the theorem implies that the divergence of the leading-twist-2 conformal 
operator with two covariant derivatives vanishes in free theory and, therefore, in full QCD 
can be expressed in terms of operators including either the gluon or the photon field 



strength. One obtains [|TT 



Tj, ^ajj, J 



20(gil 



(4.33) 



where Oj^ = y^'^-L+ D+ D+ q — |9^gcr_L+g and the explicit expression for the quark- 
antiquark-gluon operator Q^a\fjL given in Eq. (4.23) in [ll]]. Qt^i is obtained from 
Q'a\^l substitution gGi^^i, eqF^^,. The vacuum expectation value of the l.-h.s. of 

( [4.33| ) in the EM field vanishes, whereas {{Qa^^)) is equal to a constant times the quark 
condensate |jll]|. It follows that the reduced matrix element of the complicated quark- 



antiquark gluon operator Q^f^\fi i^ Held can be calculated exactly in terms of the 

condensate, and this is how Eq. ( |4.32| ) emerges. This relation is of principal importance 
since it tells that contributions of gluon and photon operators cannot be separated in a 
meaningful way. Hence the use of an approximation for gluon contributions necessarily 
entails the use of the corresponding approximation for photon contributions as well. In 
particular, the quark-antiquark-photon DAs have to be taken into account to the same 
accuracy in the conformal expansion as the quark-antiquark-gluon ones, notwithstanding 
that the former are known exactly (as given in ( [4.2U| )). 

The conformal expansion of quark-antiquark-photon DAs is constructed following the 



approach of Ref. p3. To this end one has to separate different quark-spin projections. 
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We define the auxiliary amplitudes 

(0|g(2;)7+7~^A..Mg(-^)|0) = {qq)F,,S^^iv,pz), 
{0\qizh.^+F,,ivz)qi-z)\0) = {qq)F,,S^\v,pz). 

Then 



The auxiliary amplitudes are expanded in conformal polynomials as 



(4.34) 



(4.35) 



a] 



foo + fio a 



+ flO Oig - 



a. 



+ foi{otg - 3ag 



+ foi (ttg - 3aq 



(4.36) 



As for the quark- antiquark-gluon DAs, the contributions in /oo have conformal spin-3, 
those in /oi and /lo conformal spin-4. The constants fij can be easily obtained by imposing 
the relations 



J VaS^^ia) = /oo, J VaS^^{a)ag,g = 0, 



and using the exact expressions ([4.20|) for and T^. We find 



00 



14 7 

1 ; /lO = y 5 /oi = 2 



(4.37) 



(4.38) 



With ( [4.35|) we then arrive at the following expressions: 



30a 



'1 - a„) + 



14 



ag{l 



an 



- (a? + aj) 



60al{aq + ag)(4 - 7{aq + a^)), 



30ag[aq — aq) 



7 7 

1 + - (Sag - 3) + -ag 



60al{aq - ag)(4 - 7{aq + aq)). 



(4.39) 



We would like to stress that the conformal expansion of T^" = is nontrivial: to any finite 
order in the expansion one obtains a nonzero result. This complication arises because T^, 
in contrast to the chiral-even twist-3 DAs, does not have simple transformation properties 
under the conformal transformation and corresponds to the difference of two amplitudes 
S^^ and S^^ with fixed spin projections. The symmetry S'^^ = S^^ is broken to each finite 
order in the expansion and is only restored in the infinite sum. Note that if conformal 
partial waves are taken into account, then the first A^ moments of vanish. In the next 
section we will consider an example to illustrate why using the truncated expressions in 
( [4.39| ) is mandatory for selfconsistency of the OPE. This complication was overlooked in 
pA|, O] so that the corresponding expressions have to be revised. 
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4.3.3 Equations of motion and two-particle twist-4 distribution amplitudes 



With the three-particle DAs in hand, we are in a position to calculate the two-particle 
DAs. To this end, it is convenient to use operator identities for the derivatives of the 
relevant chiral-odd nonlocal operators and take the vacuum averages in the BG field. For 
example, one obtains[3 



d _ - 
-I- — q{x)q{-x) + J dvvq{x){gG^y + eqF^y){vx)q{-x) 



dxu 



— I 



dvq{x)xp{gGp^ + eqFp^){vx)a^uq{ 



-X) 



_d_ 

dx. 



q{x)a^uq{-x) = -id^q{x)q{-x) + / dv q{x){gG^u + eqF^y){vx)q{-x 



-1 



dvvq{x){gG^^ + eqF^^){vx)a^yq{-x). 



(4.40) 



These identities can be obtained from the ones derived in ^ by the replacement 
gG gG + CqF. Taking matrix-elements and comparing with the definitions in ( |2.19|) 
and (|4.21| ), one obtains the following relations 



h^{u) 
A{u) 



— ^ / da. 
du 



duq — 

ag 



an 



(2m 



-[S{a) + S^{a)] + T,Xa) - T2(a) 



dv{Au + 1 — Qv)h^{y) 



(4.41) 



-4 



da„ 



dan 



a. 



aq — aq— {2u — \] 



a„ 



{T2{a)~T^{a))-[S{a) + S^{a)] 



from the first and from the second of the identities in ( [4.4Q|) , respectively. These relations, 
again, correspond to (A. 13) and (A. 14) in |31|] with the replacement S — >■ S + S^. Inserting 
expressions for the three-particle DAs at next-to-leading order in the conformal expansion 
and using ( [4.32|) , we obtain for h-^: 



h^{u) = -10(1 + 2h+)GI^\2u - 1). 



(4.42) 



The cancellation of contributions ~ C« is a consequence of conformal symmetry. To see 
this, notice that corresponds to the configuration with both the quark and the antiquark 
having "minus" light-cone projections and therefore a well-defined conformal expansion 
in terms of Gegenbauer-polynomials Cj^^ where j is the conformal spin, cf. 0, 0. The 
symmetry h.y{u) = h.y{l — u) implies that only odd conformal spins can contribute, and 
in addition the contribution of the lowest possible spin j = 1 vanishes by EOM. Hence 
j = 3, 5, . . ., and in particular the result in ( |4.42| ) corresponds to j = 3. The ^-contributions 

^° Here and below we use 9^ to denote the derivative of a nonlocal operator with respect to the total 
translation, cf. Il5|. 
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to three-particle DAs have j = 4; they cannot appear in the expansion of and indeed 
cancel by virtue of ( |4.32| ). 

The situation with A is different, since this function is defined off the light-cone and 
its conformal expansion is more complicated. Using the result in ([4.42|) and the second 
relation in ([4.41|) we obtain 

A{u) = AOu'^u^ {3k - + 1} + 8 - 3C2) X 

X {uu{2 + 13mm) + 2u^{10 - 15m + Gm^) Ium + 2u^{10 - 15m + Gm^) Ium) . (4.43) 

We emphasize that our results for the DAs and A are consistent with the approxima- 
tions adopted for the three-particle DAs; that is, this set of photon distribution amplitudes 
satisfies all exact QCD identities that are consequence of EOM. This issue is nontrivial. 
For example, consider an important operator identity which is an extension of Eq. (24) in 



duu 







_d_ 



q{ux)a^^q{-ux) - (/i ^ z/) V - e^^x^dt. / o^mm q{ux)^^q{-ux) 



pi pu 

-i^tiux^ jduu^ I dvq{ux) {gG^^{vx) + eqF^^{vx)}-f5q{-ux) 

Jo J~u 



+1 / ciMM / dvvqiux) {[gG^^{vx) + eqF^^{vx)]a - (n ^ i^)} q{-ux). 

Jo J-u 

(4.44) 

Taking the vacuum average in the BG field and using the definitions of photon DAs, this 
implies one more relation 



A(m) 



dv{2v - 1) h^{v) 
1 



da. 



dan 



an 



On 



an 



(2m 



an 



{T,{a) + T2{a)-T:,{a)) + S{a) 



(4.45) 



One can check that the resulting expression for A coincides with the result given in ( |4.43|) if 
and only if full expressions for the three-particle DAs and in particular the nonzero result 
for T4 (|4.39|) are used. In Ref. |]l3l the DAs and A were derived using a truncated 
version of ( 4.44| ) neglecting Tj and quark-antiquark-gluon operators. In view of relation 
( [4.32| ) this approximation is inconsistent and the resulting DAs are superseded by the ones 
obtained in this paper. 



This identity is very useful in practical calculations since it allows to express the chiral-odd operator 
with two free Lorentz-indices through a simpler nonlocal operator with one open index and three-particle 
operators of higher twist. 
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5 Summary and Conclusions 



In this paper we have given the first comprehensive analysis of photon distribution am- 
phtudes. The special interest (and complication) springs from the interplay between hard 
(perturbative) and soft (nonperturbative) contributions to the photon structure, as dis- 
cussed in Sec. 2.2. This interplay becomes nontrivial in higher twists, as QCD equations 
of motion receive additional contributions from photon operators. (Alternatively, one may 
leave equations of motion intact, but take into account contact terms.) 

We have defined and classified all two- and three-particle distribution amplitudes of 
the photon up to twist-4. For the chiral-odd ones, we could largely draw on the results 



obtained in Ref. [TT[], the remaining definitions resemble those for the p-meson discussed 
in [|3l] , ^ . The main theoretical result of our analysis is the conformal expansion of DAs 
in presence of the background field. We provide a selfconsistent approximation for photon 
DAs up to twist-4 that contains a minimal number of nonperturbative parameters and is 
consistent with the exact QCD equations of motion. The final expressions are collected 
in Eqs. (twist-3) and (^1271) , ( ^39|) , (^1421) , (|^) (twist-4). The numerical values 



of the parameters are estimated, when possible, using QCD sum rules. In particular, we 
give a new estimate of the magnetic susceptibility of the quark condensate and derive a 
sum rule for the normalization of twist-three distributions. All parameters are accessible 
to lattice-simulations, and we hope that the lattice-community will take up the task and 
provide more accurate estimates. 

We believe that the results of this work can be applied to the studies of hard exclusive 
processes involving photon emission (or absorption) in the framework of QCD factoriza- 
tion. A few examples are mentioned in the introduction. A more detailed discussion and 
concrete calculations go beyond the tasks of this paper. 
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Appendices 

A The Conformal Basis 

For a generic three-particle operator a "conformal basis" can be constructed as follows 
3^: conformal symmetry allows one to fix the total three-particle conformal spin 



J = ji + 32 + ja + of a state. We define a set of functions Yj- by requiring that, 
in addition to a fixed J, they also have a definite value of the conformal spin in the 
quark- ant iquark channel (12), for definiteness: j = ji + j2 + n with n = 0, . . . , A^. 

Taken together, these two conditions determine the polynomials Yj^'^^^ uniquely and 
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yield the following expression: 

-ii-i2 p(2i3-i,2i-i)('^ _ 2rt,l p(2ii-i.2i2-i) ( 

1- as 



Yjf\a.) = (1 - a3)— - 2a,) Pj^f "^^ ( ) . (A.l) 



Here P^'^^ (x) is the Jacobi-polynomial. 

Jj 



The basis functions Yj- (aj) are mutually orthogonal with respect to the conformal 



scalar product: 







' Va al^'-'af-\f'-' Fjf f],^^)' (a,) = Uj, 6jr , (A.2) 



where 

r(j+ji-j2)r(j-ji+j2) r(j-j+j3)r(J+j-j3) 



AT, 



J] 



r(j-ji-j2+i)r(j+ji+j2-i)(2j-i) r(j-j-j3+i)r(j+j+j3-i)(2j-i) 

(A.3) 



The above construction of the conformal basis involves the obvious ambiguity in what 
order to couple the spins of partons to the total spin J. Choosing a different two-particle 
channel one obtains a different conformal basis related to the original one through the 
matrix VL of Racah 6j-symbols of SL(2,M), e.g. 

yjf'^^i) = E ^^AJ) Yif'ix.) . (A.4) 

jl-+j2<j'<J-j3 

The properties of the Racah 6j-symbols as well as explicit expressions in terms of the 
generalized hypergeometric series 4-^3(1) are summarized in . 



B QCD Sum Rules 

Given the considerable time that has passed since the existing numerical estimates of the 
parameters of photon DAs were obtained 0, we feel that the corresponding QCD 

sum rules deserve a fresh look. Some new results are also included: we calculate the radia- 
tive correction to the leading-twist sum rule and derive a sum rule for the normalization 
of twist-3 distributions. 

B.l The leading twist-2 distribution amplitude 

The QCD sum rules for the leading-twist photon DA (f)f{u) are derived from the correlation 
function of the nonlocal tensor operator on the light-cone with the vector current. For 
massless quarks: 

1 

z J dSe-'^y{0\T[q^Mmcr,z[0,z]q{z)\0) = tz,^^{qq) J due-''^''^^u,q') (B.l) 
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with 
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6 q 
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where 6'{x) = {d/ dx)5{x) and 
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B.3) 



Here the [ ]+ prescription is defined as usual 



[/(«)]+ = /(«)-5M / dwf{w), 

Jo 

where mq = or mq = 1 is the position of the singularity, and ttiq parametrizes the mixed 
condensate: {qcrgGq) = ml{qq). The radiative correction to the quark condensate in ( |B.3|) 
is a new result]^ By expanding ( p.l|) around 2; = in terms of conformal invariant local 
operators one obtains the functions A„(g^) given in Eq. (4.7) in up to the terms in 



B.1.1 Magnetic susceptibility 

The magnetic susceptibility of the quark condensate, Xy is obtained from the correlation 
function 

) = / iu^u. = -2 + ?^ (in ^ + 1) - 1^ + . + . . . , (B.4) 



where we have also indicated that the gluon condensate contribution vanishes in factoriza- 
tion approximation pSf. We write the physical spectral density of this correlation function 
as a sum over several narrow resonances plus a smooth continuum starting at a threshold 
So- Assuming quark-hadron duality, the continuum contribution can be represented by 
the perturbative imaginary part of the radiative correction, so that 

M) - - E - 1^ («^) f (B.5) 

where fi and denote the couplings of the resonance i to the vector and tensor current, 
and are defined as in ( |3.7| ). The factor 8/3 in front of the integral is related to the 
anomalous dimensions defined in ( p.4|) : 1/2 (70 — 7(gg)) =8/3. Once these parameters are 
fixed, the magnetic susceptibility is obtained from ( p.5| ) as x = x(0)- 

In a first approximation one can retain the contribution of the lowest-lying p-meson 
state only and use the same value of the continuum threshold as obtained in the sum rules 



12 



The expression given in p3[ for the local limit z — s- contains a typo. 
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for the correlation function of vector currents, Sq 
rule estimate = 160 MeV one obtains 



1.5GeV2 m. Using the QCD sum 



X(l GeV) 



fpf, 



+ 



]a., 1 



m^{qq) 3tt Sq 



3.0 GeV" 



(B.6) 



which is an improvement over the pure VDM estimate in ( p.9|) . Including more resonances 
and at the same time increasing the continuum threshold, one may hope to get a better 
accuracy. The corresponding sum rules can be constructed in two different ways: either 
matching the expansion of ([B.5| ) for — — oo to the power series in ( plj) or perform- 
ing the Borel-transformation of both expressions and matching the results for intermediate 
values of the Borel-parameter of the order of a few GeV |]2^. Once the parameters are 
fixed, the magnetic susceptibility is obtained from ( [B.5| ) as x = x(0)- 

Following the first procedure we obtain the so-called local-duality sum rules: 



-2(gg)(/i) 
-4(gg)(/i) 



i + l^flni^-i 



37r 



4a.(/i) . rnl{ii) 
^^sl + . (G^) 

OTT 



(B.7) 



For the couplings of the ground-state resonance p(770) one can use the values quoted 
in ( |3.7|) . The three local-duahty sum rules in ( p.7| ) can be used to determine the con- 
tributions of three more resonances, so that we truncate the sums in (p.7|) after four 
terms. Experimentally, the three lowest-lying resonances p(770), p(1450) and p(1700) are 
established |26]. The magnetic susceptibility x(p) then exhibits a mild dependence on 



So and on the unknown mass of the fourth resonance m4. As an example, we show in 
Fig. x(l GeV) as function of m4 with sq fixed as Sq = {m^ + 0.1 GeV)^ (solid hne) and 
(m4 + 0.2GeV)2 (dashed line). In this calculation we use {qq){lGeV) = (-0.25 GeV)^, 



rrir 



[1 GeV) = 0.65 GeV^, cf. [gSf, and the value a,(l GeV) = 0.51. 



Neglecting for the moment the uncertainty of fif^ and the local-duality approximation 
itself, we obtain from local-duality sum rules with four resonances 



X(l GeV) = (3.08 ± 0.02) GeV" 



(B.8) 



For a typical value of m^, the breakdown of the total number in contributions from the 
different resonances and the continuum is 2.73 -|- 0.30 — 0.07 + 0.04 + 0.10, which shows 
a clear dominance of p(770). Neglecting radiative corrections altogether, we would find 
X = (3.3±0.1)GeV'2_ 

As a consistency check, we relax the condition that the known value of /i/f" be used as 
input and include only three resonances in the sums in (|B.7|) . We then find ffi^l GeV) = 
(170 ± 7) MeV for m| < sq < 5 GeV^, which is consistent with the result from Borel sum- 
rules quoted above, and x(lGeV) = (3.20 ± 0.15) GeV~^, in complete agreement with 
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Figure 3: x(l GeV) as function of m4 with sq fixed as Sq = (m4 + 0.1 GeV)^ (solid line) 
and (m4 + 0.2 GeV)^ (dashed line). 



.0125^ 

0.01 / 
0.0075 / 

0.005 ' 
.0025 

-0 . 0025 I 

1 1-5 2 2.5 3 1 i_5 2 2.5 3 

Figure 4: (a) Coupling /p'/J from the Borel sum-rule (p.9| ) as a function of the Borel- 
parameter. Sohd line: nipi = 1.4 GeV, dashed line: nip' = 2 GeV. sq = {nipi + 0.1 GeV)^. 
(b) x(l GeV) from (^) with = o with two resonances, p(770) and an effective p', as a 
function of the Borel-parameter. Identification of curves as in (a). 

Alternatively, and this constitutes the 2nd procedure, one investigates the Borel sum- 
rule 

mpfpfpC +mp,fp,fpe = 

where for simplicity we have substituted the sum over higher resonances by the contri- 
bution of one effective state p' and consider the mass mpi as a free parameter. The sum 
rule (p.9| ) can be used to extract the product /p'/J as a function of M^, sq and nipi. We 
find that for iripi > 1.4 GeV there is a reasonable stability-plateau in the Borel-parameter 
M^, as shown in Fig. |[ Varying mpi between 1.4 and 2.0 GeV and fixing the continuum 
threshold as Sq = {nip' + 0.1 GeV)^, we find from ( p.5| ) 

X(l GeV) = (3.15 ± 0.10) GeV"^ (B.IO) 

which is perfectly consistent with the value obtained from local-duality sum-rules. 

We are now in a position to quote our final value for x(l GeV). Averaging over the 
results from local-duality sum rules with three and four resonances, and the result from 
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Borel sum-rules, we obtain the central value 3.15 GeV^^. In estimating the uncertainty, 
we take into account that the p(770) is the dominant contribution, so that A^/x = 
^(/p/J)/(/p/J) good approximation. Adding in quadrature the uncertainties from 
Afp = lOMeV, A/J = lOMeV and the uncertainty quoted in (|RlO| ), we finally arrive at 
the estimate quoted in ( p.l3|) . 



B.1.2 Higher moments 



Let us now turn to the shape of the photon DA. At first sight, assuming the vector- 
dominance approximation in (|B.1|) l/(— g^) — > at 0, one is led to the conclusion 
that the photon DA schould be strongly peaked at the end-points, with the quark conden- 
sate contribution 6{u) + 6{1 — u) somewhat smoothened by contributions of condensates 
of higher dimension. We start by giving a qualitative argument why this picture is wrong. 
To this end, it is convenient to consider the Gegenbauer-moments of the correlation 



function (|B.1|) defined as 



2(2n + 3) 



3{n + l){n + 2) Jq 



2{2n + 3) 



- + 



[n + l){n + 2) ml 
6 



(B.ll) 



where for simplicity we suppress radiative corrections. Obviously Xo(q'^) = ^-^id the 

coefficients (pn in the expansion of the photon DA over the Gegenbauer-polynomials ( p.4j) 
are given by 0„ = Xniq^ = 0)/x- Assuming that the correlation function Xn for arbitrary 
n is saturated by the contribution of the p-meson and using one obtains 



iVDM 
n 



{2n + 3) 



(B.12) 



and 



= 6mm ^ i (2n + 3) C^J^{2u - 1) = - \b[u) + b{u) 



:b.i3) 



as expected. 

To see why this is wrong, consider the expansion of the VDM-type contribution in a 
power series 



The comparison of the two expansions ( [B.ll| ) and ( [B.14D shows that in order for them to 
be consistent the mass of the resonance has to grow rapidly with n: m„ ~ n. Alternatively, 
one can say that the low-mass p-meson constribution becomes irrelevant for large n. 

To see what this means, let us generalize the vector-dominance approximation by 
allowing the effective mass of the resonance to be ra-dependent and, in particular, choose 
the mass such as to reproduce the first two terms in the operator product expansion: 



,eff 



m^(n + l)(n + 2)/6, 



(B.15) 
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It is easy to see that in this approximation 



(2n + 3) ml^ff _ 2(2n + 3) 
3 3{n+l){n + 2)- 



One obtains in this model 



(B.16) 



^ 2(2n + 3) ./y, 

^ 3n + l n + 2 " ^ ' 



(B.17) 



that is a flat DA with no momentum fraction dependence: the peaks at the end-points have 
disappeared. Note that our argumentation only invokes the structure of nonperturbative 
power-like corrections in the operator product expansion. Taking into account perturbative 
gluon radiation is expected to further suppress the end-point behaviour, as illustrated in 
Fig. for the example of the instanton model. 

Trying to be somewhat more quantitative, we can write the generalization of (p.5|) to 
arbitrary Gegenbauer-moments: 



Xn{q^){qq) = - J] 



mj — 



ds 



:b.i81 



s\s — 



where (\)i^n is the nth coefficient in the Gegenbauer-expansion of the DA of the ith reso- 
nance. Note that 7„ — 7(gg) = 4Cir For n = 2, one can construct local-duality 
sum rules: 



14 



^mifiil{^)4>iM = -^{qq){^^) 
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9 vr 



6 /i 



36 



(B.19) 



where we suppress a third sum rule, analogous to the 3rd equation in ( p.7| ), because 
the factorization approximation for the gluon condensate contribution is known to be 
unreliable. 

A comparison with ( |B.7| ) shows that the expression on the r.-h.s. of the first sum rule 
is by a factor 1.36 larger (assuming sq = 2.5 GeV^), whereas the r.-h.s of the second sum 
rule increases by roughly a factor 3. This is consistent with the pattern suggested above, 
that higher resonances start to play a decisive role. Using a model with two resonances, 
p(770) and p(1450), and the continuum treshold Sq = 2.5 GeV^ one obtains the small 
value 02 = 0.07. The problem with this procedure is that the fit "wants" a negative value 
of the second Gegenbauer-coefficient 0p,2, which contradicts a more direct QCD sum rule 
evaluation of this coefficient from diagonal sum-rules ||23[. If the value p3|, 0p,2 = 0.2±0.1, 



is enforced, and, say, the mass of the resonance is taken as free parameter, the result for 
02 increases by almost an order of magnitude. In our point of view such an instability 
indicates that no reliable conclusion can be drawn on the value of 02 in the QCD sum rule 
approach. 
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B.2 Twist-3 Distributions 



To leading-order accuracy of the conformal expansion we only need the parameter f^-y 
defined in (14. 11). It can be estimated from the correlation function 



zjd'xe-''^^{0\TmG.Ml.l5qm{xhM^m = h,{q^){qz){zpq, - {qz)g,p) + . . . , 

(B.20) 

where the Lorentz-structure is chosen in order to pick up the relevant coupling of twist- 
3 corresponding to the transverse projection Gzu GzL and to suppress contributions 
of twist-4. /s^ is obtained as /37(0). The operator product expansion of fs-yiq"^), up to 
contributions of dimension 6, reads 

h^iq ) = 7ZZ^{^^'^—r + const. + ^ ( — ) - n Ta ' (B-21) 



An 7r2 \36 J 24g2 \ tt / 9 

and the p-meson contribution can be obtained by a standard QCD sum-rule calculation, 



Using So = 1.5 GeV^ and G^) = 0.012 GeV"^ 0, we obtain from this sum rule (at the 
normalization point /x = 1 GeV) 

C3p(l GeV) = 0. 10 ± 0.05, (B.23) 
which is by a factor 3 larger than the result quoted in and obtained in Ref. from 



the diagonal correlation function of two quark-antiquark-gluon operators. 

In vector-dominance approximation = —fpCsp- Including the continuum contribu- 
tions, fs-y is obtained asQ 

A,(0)^-/,\3, + ^$i (B.24) 

Using the value of Csp given in ( p.23|) we obtain the estimate (at the normalization point 
fi = 1 GeV) 

= -(0.0039 ± 0.0020) GeV^ (B.25) 

with So = 1.5 GeV^. The correction to the vector-dominance approximation is small, 
-5%. 

To next-to-leading order in the conformal spin, we have two parameters, and u}^ , 
which are defined by the following matrix-elements: 

3 
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r ~ 1 3 ~ 3 

(0|g/75 iDz, gG^^z^" q - -{idz)q^-i^gG^^z^'q\Q)F = eqi{pzfF^J;i^— + 0{zj), 

^ ^ (B.26) 



{m{9Ga.uz''{i Dz)-{iD z)gGa,,z'')q\Q)F = eq{pzfF,,h^—u'^ + 0(z,), 
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In order to obtain this formula, we express the continuum contribution to (B.21) via a dispersion 



relation ~ J^^^ ds — ji^ — sq, where /li? is the factorization scale. The quadratically divergent term ~ 



matches the contribution of low momenta in the perturbative contributions of twist- 1 in ( 2.1§| ) which have 
to be subtracted to avoid double counting. We prefer to define perturbative contributions in such a way 
as to include all momenta down to zero, and therefore have to omit quadratically divergent parts of the 
higher- twist operators. Thus we only keep the term in sq. 
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where [., .] stands for the commutator. In the vector-dominance approximation = 
and cij^ = uj^p. The numbers given in ( [4.10| ) are quoted from ||3^ and have to be taken with 
caution. Since even for leading twist we have been not able to draw reliable conclusions 
on similar parameters from the sum rules, we consider the given numbers as an order-of- 
magnitude estimate. From the general experience of similar calculations we expect that 
this is rather an upper bound. 

References 

[1] p. Kroll and M. Raulfs, Phys. Lett. B 387 (1996) 848; 

A. V. Radyushkin and R. T. Ruskov, Nucl. Phys. B 481 (1996) 625; 
I. V. Musatov and A. V. Radyushkin, Phys. Rev. D 56 (1997) 2713. 

[2] D. Miiller et a/., Fortsch. Phys. 42 (1994) 101; 
X. D. Ji, Phys. Rev. Lett. 78 (1997) 610; 
A. V. Radyushkin, Phys. Lett. B 380 (1996) 417. 

[3] M. Beneke, T. Feldmann and D. Seidel, Nucl. Phys. B 612 (2001) 25; 
S. W. Bosch and G. Buchalla, Nucl. Phys. B 621 (2002) 459; 
A. Ali and A. Y. Parkhomenko, Eur. Phys. J. C 23 (2002) 89. 

[4] E. Witten, Nucl. Phys. B 120 (1977) 189; 

W. A. Bardeen and A. J. Buras, Phys. Rev. D 20 (1979) 166 [E: D 21 (1979) 2041]; 
D. W. Duke and J. F. Owens, Phys. Rev. D 22 (1980) 2280. 

[5] J. S. Schwinger, Phys. Rev. 82 (1951) 664. 

[6] L. F. Abbott, Nucl. Phys. B 185 (1981) 189; Acta Phys. Polon. B 13 (1982) 33. 

[7] E. V. Shuryak and A. 1. Vainshtein, Nucl. Phys. B 199 (1982) 451. 

[8] M. S. Dubovikov and A. V. Smilga, Nucl. Phys. B 185 (1981) 109. 

[9] B. L. loffe and A. V. Smilga, Nucl. Phys. B 232 (1984) 109. 

[10] V. A. Novikov et al, Fortsch. Phys. 32 (1985) 585. 

[11] 1. I. Balitsky, V. M. Braun and A. V. Kolesnichenko, Sov. J. Nucl. Phys. 48 (1988) 
348 [Yad. Fiz. 48 (1988) 547]; Nucl. Phys. B 312 (1989) 509. 

[12] V. L. Chernyak and A. R. Zhitnitsky, JETP Lett. 25 (1977) 510; Sov. J. Nucl. Phys. 
31 (1980) 544 [Yad. Fiz. 31 (1980) 1053]; 

V. L. Chernyak, V. G. Serbo and A. R. Zhitnitsky, JETP Lett. 26 (1977) 594; Sov. 
J. Nucl. Phys. 31 (1980) 552 [Yad. Fiz. 31 (1980) 1069]; 

G. P. Lepage and S. J. Brodsky, Phys. Lett. B 87 (1979) 359; Phys. Rev. Lett. 43 

(1979) 545 [E: 43 (1979) 1625]; Phys. Rev. D 22 (1980) 2157; 

S. J. Brodsky, G. P. Lepage and A. A. Zaidi, Phys. Rev. D 23 (1981) 1152; 



31 



[13: 

[14 

[is: 

[16 

[17: 

[18 
[19 

[20 
[21 
[22 
[23 
[24 
[25 
[26 
[27 
[28 
[29 
[30 

[31 
[32 
[33 



A. V. Efremov and A. V. Radyushkin, Phys. Lett. B 94 (1980) 245; Theor. Math. 
Phys. 42 (1980) 97 [Teor. Mat. Fiz. 42 (1980) 147]. 

A. Ali and V. M. Braun, Phys. Lett. B 359 (1995) 223. 

V. Y. Petrov et al, Phys. Rev. D 59 (1999) 114018; 

M. Praszalowicz and A. Rostworowski, local NJL model," Phys. Rev. D 64 (2001) 
074003. 

I. I. Balitsky and V. M. Braun, Nucl. Phys. B 311 (1989) 541. 
V. M. Belyaev et al, Phys. Rev. D 51 (1995) 6177. 

J. C. Collins, Renormalization, Cambridge University Press, Cambridge (1984). 
I. I. Balitsky and V. M. Braun, Nucl. Phys. B 361 (1991) 93. 

J. D. Bjorken and S. D. Drell, Relativistic Quantum Mechanics, McGraw-Hill, New 
York (1964). 

M. Wiisthoff, Phys. Rev. D 56 (1997) 4311. 

D. Y. Ivanov and M. Wiisthoff, Eur. Phys. J. C 8 (1999) 107. 

V. L. Chernyak and A. R. Zhitnitsky, Phys. Rept. 112 (1984) 173. 

P. Ball and V. M. Braun, Phys. Rev. D 54 (1996) 2182. 

M. A. Shifman and M. I. Vysotsky, Nucl. Phys. B 186 (1981) 475. 

I. I. Balitsky and A. V. Yung, Phys. Lett. B 129 (1983) 328. 

D. E. Groom et al. [Particle Data Group], Eur. Phys. J. C 15 (2000) 1. 

V. M. Belyaev and Y. I. Kogan, Yad. Fiz. 40 (1984) 1035. 

I. I. Balitsky, A. V. Kolesnichenko and A. V. Yung, Yad. Fiz. 41 (1985) 282. 

V. M. Braun et al, |hep-ph/0206"305 . 



M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147 (1979) 385, 
448. 

P. Ball et a/., Nucl. Phys. B 529 (1998) 323. 

P. Ball and V. M. Braun, Nucl. Phys. B 543 (1999) 201. 

A. R. Zhitnitsky, I. R. Zhitnitsky and V. L. Chernyak, Sov. J. Nucl. Phys. 41 (1985) 
284 [Yad. Fiz. 41 (1985) 445]. 



[34] V. M. Braun and I. E. Filyanov, Z. Phys. C 48 (1990) 239 [Sov. J. Nucl. Phys. 52 
(1990) 126]. 

32 



[35] S. Ferrara et ai, Phys. Lett. B 38 (1972) 333. 
[36] V. M. Braun et al, Nucl. Phys. B 553 (1999) 355. 

[37] V. M. Braun, G. P. Korchemsky and A. N. Manashov, Nucl. Phys. B 603 (2001) 69. 



33 



